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We propose a conjugated filter oscillation reduction scheme for solving Burgers' equation with 
high Reynolds numbers. Computational accuracy is tested at a moderately high Reynolds number 
for which analytical solution is available. Numerical results at extremely high Reynolds numbers 
indicate that the proposed scheme is efficient, robust and reliable for shock capturing. 

PACS numbers: 02.70.-c, 05.10.-a, 47.40.Nm,52.35.Tc 

The fundamental equation for the description of complex fluid flow is the Navier-Stokes equation, for which the 
full solution is still extremely difficult in the full domain of physical interest. Burgers' equation |jj is an important 
simple model for the understanding of physical flows. Simulation of Burgers' equation is a natural first step towards 
developing methods for computations of complex flows. It appears customary to test new approaches in computational 
fluid dynamics by applying them to Burgers' equation. Jamet and Bonnerot solved Burgers' equation by using 
isoparametric rectangular space-time finite elements 0. Jain and Holla ra developed a cubic spline approach for 
coupled Burgers' equation. Varoglu and Finn |4[] proposed an isoparametric space-time finite element method for 
solving Burgers' equation, utilizing the hyperbolic differential equation associated with Burgers' equation. They 
obtain very high accuracy and numerical stability with a reasonable number of elements and time steps. Their 
method was compared with a least square weak formulation of the finite element method by Nguyen and Reynen 
|| . Caldwell et al Q further developed the finite element method to allow different sizes of the elements at each 
stage based on the feed back from the previous step. A generalized boundary element approach was proposed by 
Kakuda and Tosaka |7J . These authors tabulated their accurate results for moderate Reynolds numbers and compared 
their results with those of Varoglu and Finn jj] and of Nguyen and Reynen || . A bidimsional Tau-element method 
was developed by Ortiz and Pun for solving Burgers' equation with accurate results. Bar-Yoseph et al discussed 
a number of space-time spectral element methods for solving Burgers' equation |q|. Arina and Canuto M treated 
Burgers' equation by a self-adaptive, domain decomposition method called the x-formulation. Various finite difference 
schemes for Burgers' equation were compared by Biringen and Saati M). Recently, Wei et al [[11| have developed 
an accurate solver for Burgers' equation in one and two space dimensions. Most recently, Hon and Mao p2| have 
compared performance of their adaptive multiquadric scheme with many other computational methods. It is not our 
purpose to exhaust the literature. Despite of much effort, numerical solution of Burgers' equation is still a nontrivial 
task especially at very high Reynolds numbers where the nonlinear advection leads to shock waves. In fact, Burgers' 
inviscid shocks plague many standard computational algorithms. 

The purpose of this communication is to report a novel scheme for solving Burgers' equation for all possible values 
of Reynolds numbers. We propose a set of conjugated filters to solve the Burgers' equation. As the first and second 
derivatives are approximated by using two high-pass filters, the numerical errors of the high-pass filters at the high 
frequency region lead to oscillation near the Burgers' shock wave front. The proposed idea is to effectively eliminate 
such an oscillation by using a conjugated low-pass filter. This set of high-pass and low-pass filters are conjugated 
in the sense that they are derived from one generating function and consequently have essentially the same degree 
of regularity, smoothness, time-frequency localization, effective support and bandwidth. In the present work, all 
conjugated filters are constructed by using a discrete singular convolution (DSC) algorithm |13[, which is a potential 
approach for dynamical simulation |l4]|l5| ] and numerical computations of Hilbert transform and Radon transform. 
The DSC low-pass filter is given by 
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and its conjugated nth order high-pass filters are obtained by differentiation 
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where A is the grid spacing and a is a regularization parameter. FIG. 1 shows the frequency responses of the 
conjugated DSC low-pass filter, 1st and 2nd order high-pass filters at a = 3.2A. It is noted that all conjugated filters 
have essentially the same effective bandwidth. Wavelet multiscale analysis is utilized for adaptive oscillation control. 
A detailed theoretical analysis and justification of this scheme is accounted elsewhere Q. 



Burgers' equation is given by 

du du 1 d 2 u 

where u(x,t) is the dependent variable resembling the flow velocity and Re is the Reynolds number characterizing 
the size of viscosity. The competition between the nonlinear advection and the viscous diffusion is controlled by the 
value of Re in Burgers' equation, and thus determines the behavior of the solution. We consider Eq. (0) using the 
following initial-boundary conditions 

u(x, 0) = sin(7ra:), 

u(0,t)=u(l,t) = 0. (4) 

Cole has provided exact solution |[L7| for this problem in terms of a series expansion which is readily computable 
roughly for the parameter Re < 100. For the parameter Re = 100, the present calculations use 41 grid points in the 
interval [0,1]. The fourth order Runge-Kutta is used for the time integration with a time increment of 0.01. 

In the present DSC treatment, we use the DSC high-pass filters (||) for solving Burgers' equation p3| . The DSC 
kernel parameters are chosen as W = 35 and <r/A = 4.5 in these calculations. Both L m and L\ errors at 9 different 
times are listed in TABLE I. In an earlier work, Kakuda and Tosaka w\ tested their generalized boundary element 
method by using 100 elements, up to 6 iterations and the same time increment as ours (At = 0.01). Their results are 
also listed in TABLE I for a comparison. The errors in both methods are very small. The DSC results are from ten 
to 10 5 times more accurate than those of Kakuda and Tosaka J7]](K-T), while they were obtained using much fewer 
grid points. Note that at a late time, the accuracy of the present DSC algorithm reaches 15 significant figures. 

The numerical solution of Burgers' equation at a high Reynolds number (Re = 10 5 ) is very difficult due to the 
presence of shock j|. A direct application of the DSC algorithm using 64 grid points (N = 64) and a small time 
increment (At = 0.001) leads to a highly oscillatory results as shown in FIG. 2(a). The time integration is shown 
up to 0.5 time units and eventually collapses at a later time. The plot is given in the spatial interval of [0, 2], which 
is generated by an antisymmetric extension of the original numerical results in the spatial interval of [0,1]. The 
oscillation starts near 0.3 time units and is accumulated and amplified in late integrations. The solid line in FIG. 1 
shows the Fourier image of the result at t = 0.5. The image has two large peaks, one near the zero and another near 
the Nyquist frequency ^. 

To analyze FIG. 2(a) further, a three-scale wavelet transform is performed and the result is depicted in FIG. 2(b). 
Daubechies' biorthogonal wavelets (D7/9) jlq] are employed for the wavelet transform. At the first scale, a response 
extended over a large domain is recorded over the southeast quarter of the quadrangle which corresponds to high 
frequency oscillations in FIG. 2(a). Note that high frequency oscillations reside exclusively at the southeast quarter 
of quadrangle because the oscillations only occur at a special frequency range. The peak in the middle of the high 
frequency response region is due to the shock front, which produces similar two other narrow high frequency responses 
at the second and third scales respectively. The response of the highest amplitude at the southwest corner is enlarged 
in FIG. 2(c). Surprisingly, this part seems containing the desired solution to Burgers' equation. However, there is a 
kink in the late part of the low pass solution which does not belong to the desired solution. Obviously, had the high 
frequency oscillations been controlled in the course of integration, such a kink would not have appeared. From this 
analysis we conclude that the peak near the Nyquist frequency in the Fourier image in FIG 1 is due the undesired 
oscillations. 

The control of oscillations can be accomplished in a number of ways. For example, Godunov algorithms p9[ , 
up-wind schemes, essentially non-oscillatory (ENO) schemes [ p0|pl| ] and weighted ENO schemes p3] are standard 
methods for handling oscillations. In the present work, we propose an alternative approach which makes use of a 
conjugated low-pass filter, Eq. (Q). Since the peak of undesired oscillations is resided outside the effective bandwidth 
of the conjugated low-pass filter as shown in FIG.l, it can be removed by the filtering of the conjugated low-pass 
filter. To eliminate oscillations and preserve the true solution effectively, we design the following conjugated filter 
oscillation reduction (CFOR) scheme. We define a wavelet high-pass measure via a multiscale wavelet transform of 
the results of conjugated high-pass filter {u(xk,t n )}^ =1 (i.e. the numerical solution of Burgers' equation) at time t n 
as 

l|W"||=£||W™ ||, (5) 

m 

where || W^ || is given by a convolution with a wavelets tpmj of scale m 
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The CFOR is adaptively implemented whenever the high pass measure accesses an appropriate positive alarm threshold 
V 

||W" +1 ||-||W"||>^. (7) 

The choice of n depends on the time increment At and the grid size Ax. FIG. 2(d) shows the results under the same 
conditions as those of FIG. 2(a), obtained by using the CFOR scheme. Note that the oscillations are eliminated and 
meanwhile, a sharp shock profile is resolved. The dots in FIG. 1 shows the Fourier image of the solution at time 0.5. 
It is noted that there is little change to the image inside the effective bandwidth of conjugated filters. However, peaks 
of high frequency oscillations are effectively eliminated. 

We test this scheme for the case of Re= 10 3 and Re= 10 5 at Ax = 0.01, At = 0.001. DSC kernel parameters are 
the same as stated earlier. As shown in FIG. 3(b) and FIG. 3(c), our results are excellent. Clearly, all oscillations 
are effective removed and the shock front is very sharp. In a dramatic case, we consider inviscid Burgers' equation 
(Re=oo). As depicted in FIG. 3(d), our CFOR scheme works extremely well for this case too. The results for Re= 100 
are shown in FIG. 3(a) for a comparison. 

In conclusion, a novel approach, the conjugated filter oscillation reduction (CFOR) scheme is introduced for solving 
Burgers' equation with a wide range of Reynolds numbers. The essence of the CFOR scheme is to adaptively implement 
a conjugated low-pass filter to effectively remove the accumulated numerical errors produced by a set of high-pass 
filters. The conjugated low-pass and high-pass filters have essentially the same degree of regularity, smoothness, 
time-frequency localization, effective support and bandwidth. In this work, all conjugated filters are constructed by 
using discrete singular convolution kernels p3|| . 

The accuracy of the approach is tested at a moderately high Reynolds numbers (Re=100) for which analytical 
solution is available. While using much fewer grid points, our results are about 10 to 10 5 times more accurate than a 
previous finite element approach Ul . For extremely high Reynolds numbers (Re> 10 5 ), the DSC algorithm develops 
severe oscillations near the shock front. The CFOR scheme is proposed to minimize error accumulations and to resolve 
the shock front. It is found that the present scheme is very accurate and robust for integrating Burgers' equation over 
all possible Reynolds numbers. 

We note that the present approach is very general. It can be applied to the numerical solution of other partial 
differential equations, particularly, compressible flows and hyperbolic conservation laws. Moreover, the CFOR scheme 
can be implemented along with any other standard computational methods, such as high-order central difference 
schemes, finite clement methods and spectral approximations. A theoretical analysis of the present approach and an 
evaluation of various wavelets and filter banks for the CFOR scheme will be reported elsewhere [[16| . 

This work was supported in part by the National University of Singapore. The authors thank Professor C.-W. Shu 
for useful discussions about the ENO and WENO schemes. 
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TABLE I. Comparison of errors for solving Burgers' equation 





K-T 


DSC 


Time 


Loo 


^oo 


Li 


0.4 


2.6(-02) 


2.4(-03) 


2.2(-04) 


0.8 


2.9(-02) 


3.3(-03) 


2. 9 (-04) 


1.2 


1.8(-02) 


4.7(-04) 


1.1 (-05) 


3.0 


6.9(-03) 


7.6(-08) 


1.1 (-08) 


10.0 




3.1(-11) 


1.2(-11) 


30.0 




1.4(-12) 


8.6(-13) 


60.0 




6.9(-14) 


4.4(-14) 


90.0 




3.6(-15) 


2.3(-15) 



Figure Captions 

FIG. 1. Graph of the frequency responses of the conjugated DSC filters (in the unit of 7r/A). The maximum 
amplitude is normalized to the unit. Stars: conjugated low-pass filter; Dashed line: 1st order high-pass filter; Dash- 
dots: 2nd order high-pass filter. Solid line: Fourier image of the numerical solution of Burgers' equation (t = 0.5, 
Rc= 10 5 ) with oscillations; Dots: Fourier image of the numerical solution obtained by using the CFOR scheme. 

FIG. 2. (a) The oscillatory numerical solution of Burgers' equation (Re= 10 5 , At — 0.001,7V = 64, t — ~ 0.5); 
(b) Three scale wavelet analysis of FIG. 2(a): (c) The last scale low pass response of FIG. 2(b); (d) The CFOR 
solution of Burgers' equation (Re= 10 5 , At = 0.001, TV = 64, t = - 0.5). 

FIG. 3. The CFOR solutions of Burger's equation at t =0.2 (i), 0.5 (ii), 1.0 (iii), 1.5 (iv) and 2.0 (v) (At = 
0.001, N= 101). (a) Rc= 100; (b) Re= 10 3 ; (c) Rc= 10 5 ; (d) Rc= oo. 
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